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Abstract 

We study the decrease of fluctuations of diagonal matrix elements of observ- 
ables and of Husimi densities of quantum mechanical wave functions around 
their mean value upon approaching the semi-classical regime (h — > 0). The 
model studied is a spin (SU(2)) one in a classically strongly chaotic regime. 
We show that the fluctuations are Gaussian distributed, with a width a 2 
decreasing as the square root of Planck's constant. This is consistent with 
Random Matrix Theory (RMT) predictions, and previous studies on these 
fluctuations |], ||. We further study the width of the probability distribu- 
tion of fi-dependent fluctuations and compare it to the Gaussian Orthogonal 
Ensemble (GOE) of RMT. 



1 e-mail philippe.jacquod@iph.unine.ch 
2 e-mail amiet@iph.unine.ch 



1 



The behaviour of quantum mechanical wave functions in the semiclassical 
limit has recently attracted much interest. It is motivated by the fact that 
the spectrum alone cannot contain the whole information on the system. 
Roughly, one can say that in integrable systems the eigenfunctions condense 
on classically invariant torii, while in chaotic ones, where such classical struc- 
tures have been destroyed, they tend to spread uniformly over the whole clas- 
sically allowed region. Few analytical results have been obtained however in 
chaotic regimes, the most important of which perhaps is the Shnirelman the- 
orem. One formulation of this theorem would be that in the limit h — > 0, 
almost all the diagonal matrix elements of almost all quantum mechanical 
observables converge weakly to a constant over the classically chaotic region 
[f|]. A few years ago Feingold and Peres and more recently, Eckhardt 
et. al. IJ have studied the rate of this convergence for autonomous sys- 
tems where the semiclassical limit is, according to the Shnirelman theorem, 
the microcanonical phase-space (i.e. classical) average. As they mentioned, 
the "almost all quantum mechanical observables" in this formulation exclude 
projection operators, and in general all operators without smooth classical 
limit. Moreover, the "almost all diagonal matrix elements" still leave room 
for scarring of eigenstates by short periodic orbits ||. F° r those states, the 
limit can be dramatically different from the Shnirelman-predicted one. Their 
conclusion is that in a strongly chaotic system and for a smooth classical ob- 
servable A(p, q) with which a quantum operator A, Ajk '■=< j\A\k >, can be 
associated, the fluctuations of the diagonal matrix elements 



have the same order of magnitude as the mean square of the off-diagonal 
terms < |^4jA;| 2 > and decrease proportionally to the inverse of the Heisen- 
berg time 1/Tjj ~ h upon approaching the semiclassical limit, in agreement 
with Random Matrix Theory (RMT) predictions. They related the propor- 
tionality coefficient to the autocorrelation function of the classical dynamical 
variable A, C(t) := lim^oo ^ J T A(t + t)A(t)cIt, i.e. 



< Ff >:=< {A n - {A}) 2 > 



(1) 



around the semiclassical microcanonical average 



{A} = J A(p, q)5(E - H(p, q))d d pd d q/ J 5(E - H(p, q))d d pd d q 



(2) 




(3) 
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In particular, almost all diagonal elements Ajj tend to the semiclassical mi- 
crocanonical average as % — > 0. Eq.(3) states among others that quantum 
fluctuations are proportional to classical correlations. Their argument goes 
as follows : According to Shnirelman's theorem, the diagonal matrix element 

< Ej\A(t)A(0)\Ej >^ C(t) h^O (4) 

On the other hand, this matrix element is 

< EjlA^AWlEj > = ]Texp [i(E j - E k )t/h] \A jk \ 2 

k 

= 5>xp -E k )t/h] |A, fc | 2 + \A j3 \ 2 (5) 

Thus we have 

J2 exp [i(Ej - E k )t/h] \A jk \ 2 - C(t) - {A} 2 h^O (6) 

Defining the Fourier Transform of the autocorrelation function S(u) : = 
X^o C(t) exp(—iu>t)dt we have 

C(t) - {A} 2 } dt E 3 - E k (7) 

-oo L J 

Then, under the assumption that as Ej —>■ E k , the eigenfunctions \Ej >, 
\E k > and |± >:= ^{\Ej > ±\E k >) are qualitatively similar, i.e. : 

A jk ^< -\A\ + > 
{A} ~< +|y4|+ >~< — 1^4| — > (8) 

we have 

A jk ^< -\A\+ >= ^{Ajj - A kk + A jk - A kj ) (9) 

Finally defining the fluctuations as Ej := Ajj — {A} and assuming statistical 
independence of the F/s, i.e.: < F 2 >=< E 2 >= 2 < \A jk \ 2 > we get 
eq. (3). Illustrations of this result on the double rotator model ||, the 
bakers map and the hydrogen atom in a magnetic field fl]] nicely confirmed 
these predictions. These are to our knowledge the only works that delt 
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with the qualitative description of the approach to ergodicity of quantum 
mechanical wave functions. Here, we extend these results to a kicked (e.g. 
non autonomous) system. We will focus on the fluctuations of the Husimi 
density of the eigenstates, i.e. study the fluctuations of the diagonal matrix 
elements of the projection operator over coherent states ||. The Hamiltonian 

H:=^-Sl + ^S y g 5{t-nT) (10) 

is expressed in term of the SU(2) spin operators S x , S y and S z , while < 
k < 2tc. Models of this kind have been extensively studied || and are usually 
referred to as "kicked tops" . They alternatively represent a spin which evolves 
during a time T under the influence of an integrable hamiltonian after which 
it undergoes a rotation of angle k around the y-axis. It thus defines the time 
evolution (Floquet) operator : 

U T :=eM-i^S y )eM-^S 2 z ) (11) 

The above argument leading to eq. (3) must be slightly modified in order 
to apply to the map defined by eq.(10) and eq.(ll). Instead of working 
with energy eigenstates \Ej > of an autonomous Hamiltonian, we deal with 
quasienergy eigenstates \ojj > of an unitary time evolution operator. As 
a consequence, the microcanonical average of eq.(2) is replaced by a phase 
space integral restricted to the corresponding connected chaotic region. In 
our case and in a strongly chaotic regime eq.(2) reads : 

{A} = I A(e, <p) sm(0)d6d<f>/ I sm(9)d9d(p = — [ A(9, <p) sin(0)d0d<#12) 
Js 2 Js 2 An Js 2 

i.e. we integrate over the whole sphere S 2 instead of the energy surface. In 
the semiclassical limit, the diagonal matrix elements 

< uo^A^A^ujj >= J^exp - w k )t/h] \A jk \ 2 -> C{t) (13) 

k 

provided the regime studied is classically strongly chaotic. Moreover, a sim- 
ilar argument as before leads to 

lA^nSd^-u^/h) (14) 
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and hence we recover eq.(3). Here, we concentrate on the study of the eigen- 
states of the unitary operator eq.(ll) in the regime T = 50 and k = 1.2. 
By standard numerical computation of the Lyapounov exponent |9| over the 
whole phase space, we checked that in this regime the classical motion is 
strongly chaotic. Moreover, we checked that the quantum mechanical oper- 
ator eq.(ll) exhibits the usual characteristics of quantum chaos : its level 
spacings statistics and spectral rigidity follow the predictions of the Circular 
Orthogonal Ensemble (COE) of RMT. 

As mentioned in [jl| the Shnirelman theorem leaves room for wave functions to 
show large deviation from the semiclassical limit value. It only states that the 
proportion of such wave functions should be negligible, i.e. in the semiclassi- 
cal limit, they build a subset of zero measure. For "almost all eigenf unctions" 
then, the variance of these fluctuations should vanish as % — > 0. However, 
this decay can be substantially perturbed by scarring of eigenfunctions by 
short periodic orbit [3| : Scarred eigenfunctions are front-line candidates for 
exceptions to the Shnirelman theorem ! Thus they could significantly - and 
negatively - affect our results. It has been suggested that scarring manifests 
itself in deviations of RMT predictions in the level curvature distribution [|7| 
[§]. Though not yet rigorously proven, this statement is now widely accepted. 
This distribution for the model defined by eq.(ll) in the regime studied is 
shown in fig.l. The remarkable agreement with RMT (full curve) prediction 
indicates a small number of scarred eigenstates, an agreement which was al- 
ready obtained on a similar model 0. Therefore, scarring is not likely to 
influence our study. 

Let us briefly outline our method. Our goal is to study the behaviour of 
eigenstates of eq. (11) in the semiclassical limit, i.e. as h — > 0, S — > oo so 
as to leave the product hS constant. A peculiarity of such systems is that 
the parameter governing the convergence to the semiclassical limit governs 
too the number of states 2S + 1 ~ 1/h and the density of states. In order to 
determine the implication of this peculiarity on our study, we will therefore 
check the validity of our results on GOE matrices. 

The Husimi density of an eigenstate \tu >= J2n=-s ^/J/- 4 > °f ec f(H) is 
defined as the projection of this state onto a coherent state \9,<p > of the 
spin SU(2) group : 

:=|<^|M>| 2 
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\0,4> >■= J(,%) sm(l) s -» cos^y+^-^lfi > (15) 

The are smooth functions of 9 and <fi and thus can be expanded in term 
of the spherical harmonics : 



0) = £*,m J^-MfYUO, <f>) (16) 



21+1 hm 1 l,m 

where I = 0, 1, 2, ...2S and m = —I, — l+l, —1 + 2, ...I. We used the convention 
to introduce the square root in this expansion. This multipole expansion 
allows us to interpret the Qf m in term of magnitude of fluctuations of size 
~ — tt in the ^-direction and ~ -Mr in the 6 direction. We will thus get 

m+l t l+x o 

quantitative results on the decrease of fluctuations as a function of their size. 
Let us recall that the Shnirelman theorem implies that as H = 1/S — > 0, 
fluctuations of fixed and non zero I must vanish. However it does not say 
anything about the behaviour of, say, £lftg\ m (g) as S — > oo when l(S) and 
m(S) are monotonously increasing functions of S , i.e. investigating such 
moments could lead us to different conclusions than that of |l|, p] . 
Using the resolution of unity : 



2s 



J d6d<f) sin 6\6, >< 9,<f>\ (17) 



47T 

the normalization condition reads : 

1 =< u\l\u >= (2S + l)£l$ fi 

=* °& = M^TT) < 18) 

i.e. the t h moment decrease as 1/S ~ % on approaching the semiclassical 
limit. In the following, we therefore divide all higher moments f2f m by f2go 

to consistantly study their decrease and introduce the notation Clf m := jjg 1 . 

Let us note that this 1/S behaviour of the Shnirelman limit is a consequence 
of the overcompleteness of the coherent states representation. On the other 
hand we have : 

cif jm = 4(25 +i) e <w-ir^ iQ c£i,_ M , m (19) 

H=— s 
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which gives a check of our numerical computation for small S. However, the 
numerical difficulty for computing the Clebsch-Gordan coefficients Cu+^.-tt.m 
for large S leads us to use the following numerically more stable and faster 
method to compute the moments Clf m . We define : 



Mj? m (uj) := 4(25 + 1) < uj\S k z S m \uj > / S k+n 



(20) 



It is straightforward to see that there is a linear relation between the M^ r 
and the Qf m : (We use the shorter notation 7 = e*^tan(6 l /2)) 



Tr 

Qk-\-m 

25 + 1 

4irS k+m 
25 + 1 

AixS k+m 



uo >< u\s K z s r : 

J d 2 7 < 7 |w >< u\S k S^ > 



(21) 



The curved letters S stand for classical quantities, and for any function /( 7 ) 
we have defined the product fTPf : 

/( 7 ) o S z := {S z - 7^)/(T) 

/(7)o5_:=(5_ + ^)/( 7 ) 

This allows us to write < j\uj >< a; 1 5^5™ I7 > as a differential operator 
acting on fi®( 7 ). The trick is then to partially integrate this expression. 
After a little algebra we reach : 



25 + 1 



47rS k+r 
(25+1 + m)! 

(25)!2 fc + m 47r Jo 



d 2 7 0f (7) o {S z o)\S^o)' 



J W d<p J due-^(l - u 2 Tl 2 Vl m {u)h* (22) 



where 



= ((25 + 2 + m)u - m - (1 - u 2 )-) k l := £ p^f 

l =771 



(23) 
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in term of the Legendre polynomials P™{u) and u = cos(0). We finally get : 

(u) - (gg + l + ^ y 1 nS f24) 

^,m^J - ( 2> S)!(25) fc + m ^21 + 1 l > mPk [ ] 

It is thus possible to obtain the Clf m through a matrix multiplication of the 
moments M§ m {uj) : 

M s {u) = MQ S (25) 
where we defined {M s {u)) Km = M^ m (u), (Q s \ m = (2S ffr )! ffi TO and 

(•A4)fc,Z = (2S)!(2S) fc 2l+T^fe- 

Numerical invertion of this last matrix allows us to get the multipoles Ctf m 
from the numerical computation of the moments Mj? m (uj). The advantage 
of this method against the direct computation of Husimi densities is the 
numerical stability. Moreover, if we are interested in the few first moments, 
say up to I <C S, then only the diagonal matrix elements up to Mf m are 
necessary. 

Fig. 2 shows a plot of a moment distribution P(Qf^) obtained through com- 
putation of 2404 diagonal matrix elements from 4 unitary matrices defined 
by eq.(ll) [] close to the regime T = 50 and k = 1.2. The agreement with the 
Gaussian fitting is remarkable and allows us to conclude that the fluctuations 
of the Qf m obey the probability distribution 

P(fl? >m ) oc exp(-(n£ m - figJ7(2o*J) (26) 

where the mean value Clf° m is the Shnirelman limit. This distribution narrows 
itself as h — > 0, until finally the "almost all" wave functions, i.e. those who 
obey the Shnirelman theorem, have converged to their semiclassical shape 
&Tm = 0, / 7^ and f2^° = 1. In other words, of m decreases as S increases. 
This decay follows a power law as shown in Fig.3. We have : 

a\ m ~ S- 1 ' 2 V/ ^ (27) 

As already mentioned, this law is valid for fixed I and m in the regime I, 
m <^ S. 

3 We have considered only the projection of (2) on even states, i.e. states which are left 
invariant by the parity >= | — /a >. 
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We further did the same study on GOE matrices. The result is shown in 
Fig. 4 and indicates a decay of the width of the Gaussian distribution of 
fluctuations of the Husimi density of the form (25). Let us note at this stage 
that the relationship between this width and the fluctuations studied in [p], Q 
is : 

(aU 2 ~< Ff > (28) 

We thus get the same 1/S decay of the fluctuations. As for the shape of 
these fluctuations, the diversity of models studied up to date leads us to 
postulate that quantum mechanical systems with strongly chaotic classical 
counterpart have gaussian distributed fluctuations of their diagonal matrix 
elements around their microcanonical classical average (Eq.(2) or eq.(12)). 
Apparently, the width of this gaussian decays like h as h — > 0. This postulate 
is to be taken with the "almost all" Shnirelman restrictions and excludes of 



course models like the kicked rotator [0 , where quantum interference effects 



lead to localization of the wave function, thus destroying the ergodicity of 
the quantum wave function []. 

Up to now we have shown that our model matches in every respect all the 
features of a GOE random matrix : its spectrum exhibits level repulsion its 
level curvature statistics correspond to the RMT predicted distribution and 
the statistical distribution of the components of its eigenvectors tends to the 
semiclassical average in the same way, which in its turn implies a decay of 
the width of the Gaussian distribution of the moments Clf m defined in (16). 
However as has already been said, there is absolutely no reason to expect a 
similar decrease of such moments when / is not small compared to S. We 
therefore turn our attention to the behaviour of such moments. 
We concentrate on the questions : 

• Is there a similar power-law decay for Clf, S \ m( - s % when l(S) and m(S) 
are increasing functions of S ? 

• Are there possibly restrictions on l(S) and m(S) for this power-law to 
stay valid ? 

Answering this questions gives us information on the minimal size A^ m of 
the relevant fluctuations. From the Heisenberg uncertainty principle we have 
a lower bound for the fluctuations size A^ m = jp^S—^ ^ h ~ 1/S and 

4 However restriction of quantum means to phase space region smaller than the local- 
ization length should lead to a similar behaviour. 
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thus an upper bound for I and m : l,m < S. For the sake of simplicity 
we will restrict ourselves to the study of m = moments, and will study 
moments with I ~ S and Vo using formula (19) with random eigenfunction 
components uj^ which corresponds to the GOE case []. 

We show the result of this study on fig.5 for l(S) = 5/2, 35/4, S and 5*9/4. 
Obviously, these S'-dependent moments decay faster than those with fixed 
I and m. Moreover a S c is likely to exist for each l(S) above which the 
magnitude of the corresponding fluctuation decays faster than a power law, 
possibly exponentially. However this latter conclusion is to be taken carefully 
because of the restricted 5-range of fig.5 Q On the other hand the I = VS- 
moment decay as a power law ~ S~ 3 / 2 , at least in the studied range of 
variation of S. 

In view of this, we conclude that, in the GOE case, the critical value l c below 
which the fluctuations are relevant either tends to a constant, or to infinity 
slower than S, i.e. 

l c ~ S a < a < 1 (29) 

On the other hand, previous study of the kicked top emphasized the quasifrac- 
tal structure of the Husimi density of its eigenfunctions in the chaotic regime 
[O]. This means that fluctuations in both directions of phase space are 
present up to the smallest scale allowed by the Heisenberg uncertainty, i.e. 
up to a size O^h 1 ^ 2 ), which is consistant with eq. (29) with a = 0.5. The fact 
that the moment Q also shown on fig.5 decays more or less as a power-law 

- VS- 3/2 (30) 

corroborates this conclusion : moments up to I ~ are relevant, i.e. a = 
1/2. 

Nevertheless, nothing forces the eigenstates of a quantum chaotical model to 
match those of a GOE matrix up to the smallest scales. It would therefore be 
highly desirable to get a condition on a like eq.(29) for a quantum chaotical 
system. This could be achieved by direct computation of Q^/s m ( S ) using 

5 The w M 's are random up to the normalisation condition 5Zu=-s l^pl 2 = 1 an< l the 
IT-parity : w M ^ either for \i = -S, -S + 2, -S + 4, ... S or [i = -S + 1, -S + 3, ... 
5-1. 

6 This restriction is due to the computation of the Clebsch-Gordan coefficients. 
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eq.(19). However, the numerical difficulty associated with the computation 
of high-order Clebsch-Gordan coefficients renders this task hardly fulfillable, 
as can be seen on fig.6 where we show results obtained for through 
eq.(19) averaged over more than 40000 states for each point. On one hand, 
the semiclassical randomness of the eigenstates is not attained for small S, 
while on the other hand, the Clebsch-Gordan coefficients limit the maximal 
spin magnitude. In other words these two effects dramatically affects fig.6 left 
and right. Considering the size of our statistics, we attribute to these effects 
the somehow erratic behaviour of ^^q- On fig-6, the solid line indicating 
a S — 1.5-decay is shown as eye-guide, and constitutes in no way a serious 
result. 

In conclusion our study of the Husimi density of eigenstates of the quan- 
tum spin system defined by (10) and (11) has confirmed the gaussian shape 
of fluctuations around the semiclassical limit. Moreover, the decay of these 
fluctuations follow the same power law as in previous studies [p], 0, indicating 
perhaps universality. While GOE results tend to confirm the quasifractality 
proposed in [|13[], numerical difficulties forbade us to check it for the quan- 
tum dynamical system. Investigations to overcome this difficulty are on their 
way. For the time being, let us just point out that the fact that GOE eigen- 
states seem to exhibit this quasifractality renders it a direct consequence of 
the randomness of the states. The maximal randomness is then bounded by 
Heisenberg's uncertainty, but beside that, the quasifractality of the states 
seem to contain no physical content. 
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part by the Fonds National Suisse de la Recherche Scientifique. 
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Figure Captions 



Fig.l: Distribution of level curvatures for the eigenstates of (11) and S = 
200, T = 50 and k = 1.2. From the remarkable agreement with RMT predic- 
tions we conclude that the ratio of scarred eigenfunctions is very small (see 
[0, §]), and should therefore not influence our study. 

Fig. 2: Moment distribution P(Clf ) as defined in (16) for a spin S = 600. 
The statistics has been computed from 2404 even states of four realisations 
of (11) taken around T = 50. and k — 1.2. The agreement with a gaussian 
(solid line) is remarkable. On inset we show the same curve on a semi-log plot. 

Fig. 3: Log-log plot of the width of the gaussiann distribution of moments 
P(Clf m ) for model (11), m = and I = 1 (squares), I = 3 (diamonds) and 
/ = 5 (triangles) vs. the magnitude of spin S. Inset shows the width of 
P(Re(f2f m )) for m = 2 and / = 2 (circles), I = 3 (squares), I = 4 (diamonds) 
and / = 5 (triangles). In both cases, the solid line indicates the S~ l l 2 decay. 

Fig. 4: Log-log plot of the width of the gaussian distribution of moments 
P(Clf m ) for GOE, m = and / = 1 (circles), / = 3 (squares) and I = 5 (dia- 
monds) vs. the magnitude of spin S. The solid line indicates the S* -1 / 2 decay. 

Fig. 5: Log-log plot of the width of the gaussian distribution of moments 
P(nf ttn ) for GOE, m = and I = S/2 (circles), I = 35/4 (squares), I = S 
(diamonds), I = 55/4 (triangles) and I = (empty diamonds) vs. the 
magnitude of spin S. The upper and lower solid lines indicate a decay of 
S^ 1 - 5 and S 1-36 respectively. 

Fig. 6: Log-log plot of the width of the gaussian distribution of moments 
p (&i(s),o) for model ( n )> l (S) = y/S and T = 50. and k = 1.2, m = and 
S = vS vs. the magnitude of spin S. The solid line indicate a decay of S 1 " 1 / 2 . 
We attribute the rather erratic behaviour of the datas to the numerically in- 
stable computation of high-order Clebsch-Gordan coefficients (see text). 
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